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Abstract 

We present simple analytic expressions for the distributions of the Bhabha scatter- 
ing process with emission of one hard photon, including weak boson exchanges, and with 
longitudinal polarization of the initial electron. The results from the Monte Carlo gen- 
erator BHAGEN-1PH, based on these expressions, are presented and compared, for the 
unpolarized case, with those existing in literature. 
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1. Introduction. 

The cross section for e + e~ i— > e + e~7 in QED has been calculated for different pur- 
poses and various distributions have been produced [1-16]. The availability of high energies 
(TRISTAN, LEP, SLC) and longitudinally polarized electron beam (SLC) makes it neces- 
sary to include weak-boson and polarization effects. For the unpolarized case a relatively 
simple expression for the square of the matrix element was obtained [2,3], for the polar- 
ized case it was advocated [4] that it is easier to use just a matrix element and not its 
square. This is true for the general case with multiphoton emission, however for the ex- 
perimentally interesting case of longitudinally polarized electron beam, by careful choice 
of the variables, in [5] was obtained a reasonably compact expression for the square of the 
matrix element. In this work that expression is improved with some other terms relevant 
on the Z boson peak and with all the relevant mass-corrections for the configurations in 
which the final fermions have angles with the initial direction larger than, say, 1 mrad. For 
extremely forward final fermions the mass-corrections reported here are not sufficient as 
discussed in [1,5,6,7,8]. We present the complete analytical form of the differential cross- 
section, omitting the details of the calculation; the expression is simple enough and can be 
easily cross-checked. Because of the experimental interest we also implemented it within a 
Monte Carlo program, allowing for a fast event generation of the process. Furthermore the 
procedure do not require the introduction of new peaks as frequently occurring in spinor 
techniques. Some results and the numerical tests are given in Section 3, while the details 
of the structure of the program are presented elsewhere [9] . 

2. The distributions. 

We consider the process 

(2.1) e+(p+) +e"(p_) ^e + (q+) +e"(g_) +-y(k). 

We define 

^ s = (P++P-f, t = (p + -q+) 2 , u = (p+ -q-f , si = (q+ + q-) 2 , 

ti = (p_ - q-f , u x = (p_ - q+) 2 , k± =p±.k , h± = q±.k, 

and 

(2.3) V = e^paP+P^q+ql , e i 2 3 = 1 • 

In the following a is the fine structure constant, Pl is the degree of longitudinal polarization 
of the initial electron (Pl = 1 for pure right-handed states) and Cy,C A the Z couplings 
to the electron. In the standard electroweak theory one has in particular 

c v = - A 1 r 4 / n29 ^ , c A = 1 



4 sin $w cos $w 4 sin 9w cos &w 

where 6\y is the weak mixing angle. We define, for later convenience, also the combinations 
(the factor 2 in C\ z was absent in [5], due to a misprint) 

^ 2 5 ^ C jz = Cy ±C 2 A , Czz = (Cy + C 2 A ) 2 ± AC V C 2 A , 

C® z = CyC A , C° zz = 2CyC A (C v + C 2 A ) . 
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With M z and T z the mass and width of the Z, 6 the usual step function and 



(2.6) 



D{s) = (s - M z f + M z T 2 z Q(s) , 



we define the following functions (slightly different from [5], to keep the expressions as 
compact as possible) 



F+(s,t,a) = 1 + 



(2.7) 



S (s - Ml) + t(t - Ml) 



+ 



D(s) D(t) 
8t [(s - Mp(t - Ml) + M|r|8( 8 )8(t)] 
D(s)D(t) 



(c+ z -2P L C° z (l + a)) 

[C+ z -2P L C zz (l + a)) , 



F~(s,t) = 1 + 



(2.8) 



s\s 



Ml) 



+ 



t(t 



Ml) 



C, 



+ 



D(s) D(t) 
st [(s - M%)(t - M 2 Z ) + Mf r|e(s)9(£)] 
D(s)D(t) 



c 



ZZ 5 



E+(s,t) =M Z T; 



(2.9) 



+ 



st (0(s)(t 



sS(s) 
Ml) 



D(t)J 
-Q(t)(s 



Ml)) 



D(s)D(t) 



1 z 



C°zz ~ ^PlCzz 



(2.10) 



E-(s,t) = -^M z T z 



sS(s) 



(t) J °^ 



+ 



D(s) D 
st (e(s)(t-Ml)-e(t)(s-Mp) 
D(s)D(t) 



C 



ZZ 



The differential cross section for the process (2.1) can be written as 



(2.11) 



da 



a" 



2tv 2 s 



(X + Y + Z) 



d 3 q + d 3 q- d 3 k 
E-\_ E_ E~, 



S 4 (p + +p- - q + - q- - k) 



where E + , E_, E 1 are the energies of the final positron, electron and photon respectively. 
The quantities X, Y, Z refer to the annihilation, the coulomb and the interference part of 
the square amplitude respectively and, omitting the terms which are numerically irrelevant 
at the interested energies and for the kinematical configurations considered (final fermions 
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not too forward), can be expressed in the form 
(2.12) 

X = [F~(s, Sl )(t 2 + t 2 ) + F+(s, s u 0)(u 2 + u 2 )] / 



Assa 



U Ui 

+ 



tl 



k+h- k-h+ k+h+ k-h. 



+ [F~( Sl , Sl )(t 2 + t 2 ) + F+( Sl , Sl , 0)(u 2 + ul)] 



Asik+k 



+ [F-(s,s)(t 2 + tj) +F+(s,s,0)(u 2 + ul)} 



1 



Ash+h. 



+V [E-(s, Sl )(t 2 - tl) + E+(s, Sl )(u 2 - ul)] 



s - Si 



2ssik+k-h+h. 



2s 2 

2 

m e 

'2s 2 



F-(s,s) 
F-( Sl , Sl ) 



t 



(h+) 2 {h.y 



(h+) 2 (h-) 2 



t \ I , . til j , . tl 



(2.13) 

y = [F~(t, ti)(s 2 + s?) + F+(t, h, 0)(u 2 + ul)] 
- [F~(t, t)(s 2 + sl) + F+(t, t, 0)(u 2 + ul)] 



1 



Att x 
1 



U Ui S Si 



k + h- k-h + k + k- h+h- 



Atk-h. 



[F-(t 1 ,t 1 )(s 2 + sl) + F+(ti, ti, 0)(u 2 + u 2 )] 



1 



Atik + h 



m e 
2lf 

2t 2 



F~(ti,ti) 



+ 7 ^)+F+(ti,ti,0) 



u 



+ 



ui 



u 



(k- 



(2.14) 



Z = 



u 2 + ul 



F+(s,t,0) 



st 



u s 

+ 



k-h + h+h- k-h. 



+ 



F+{s,ti,0) f ui 



sh 



k+h- h+h- k+h 



F+(si,t,0) ( ui si 



Sit 



k+h- k+k- k-h 



t \ F+(si,ti,0) 



Siti 



+ii + 
ti 



k-h+ k+k- k+h. 



u Si 

+ 



+V 



E+(s,t)k+ E+(s,ti)k_ E + (s u t)h+ E+{si,ti)h. 



st 



+ 



sti 



Sit 



Siti 



U — "UT 



k+k-h+h- 



—±F+(s,t u 0) 
sti {h+y siti 



9 9 9 

F + (s U tl,0)—±r 



(k+y 



m 2 . ^ + , x ul TYl 2 I . . u 2 



where m e is the electron mass and a = 8 * • The mass-terms (i.e. the terms proportional 
to ?rig) included here are the relevant ones when the final fermions have angles with the 
initial direction larger than, say, 1 mrad. For even smaller angles a few more terms have 
to be included as discussed in [1,5,6,7,8]. The contributions to Eq.s (2.12-14) proportional 
to the pseudoscalar of Eq. (2.3), not reported in [5], are relevant at the Z boson peak. 
The unpolarized cross section can be easily recovered by putting Pl = 0. 

We have verified that the unpolarized annihilation part agrees with the result of 
Ref.[17] for e + e~ \— > / u + / u~7, in the limit of massless fermions. The coulomb part can 
be verified by noting that it must be reproduced by crossing from the annihilation, i.e. 
by the exchange of with — q_. Finally, after some algebra, one can recover from our 
expressions the QED result of Ref.[2] in the appropriate limit; as a consequence of the 
presence of the ^-exchange terms, our formulae do not posses the factorization observed 
there. Such a factorization is recovered only in the unrealistic Mz = limit, suggesting 
that the connection between hard and soft photon emission observed there is related to 
the absence of a mass scale. Our results reproduce also the results of Ref. [3] for the 
unpolarized case (in the t-channel Z boson contribution we put to zero the absorptive 
part, i.e. the part proportional to the Z width). 

3. Numerical results and comparisons. 

On the basis of the above formulae we have constructed a Monte Carlo event generator, 
BHAGEN-1PH. The algorithms used in the generation and a long write-up can be found 
in Ref. [9]. However an indispensable part of the construction of the Monte Carlo event 
generator is its test, which assure that the algorithms used for the generation are well 
constructed and applied. Numerical tests check also a number of bugs in the program. The 
tests through comparisons with existing independent programs where possible only for the 
situation where the beams are unpolarized as no program which can provide cross section 
for polarized beam(s) is available to our knowledge. However as the algorithm used in 
BHAGEN-1PH for generation in the situation with the polarized electron beam is identical 
to the one used for unpolarized beams, no further tests of the algorithm are necessary. On 
the occasion of the estimation of the theoretical error of the luminosity measurement at 
LEP [10,11], extensive tests of the QED only t-channel part of the program were done, as 
BHAGEN-1PH was included as a part of BHAGEN95, a Monte Carlo program to calculate 
the cross-section for Bhabha scattering, completed with radiative corrections. There the 
accuracy of 3 x 10 -4 was reached in the 0(a) comparisons; no deviation up to this accuracy 
was observed [10,11] from OLDBIS [12], SABSPV [13], NNLBHA [14], BHLUMI [15]. To 
complete this tests we present here a comparison with OLDBIS, where only contributions 
to the cross section of the hard photon spectrum are separately compared in Table 1. This 
specific tests were not done in [10,11] as only the total error on the complete Bhabha 
scattering was relevant there. The event selection used here is the BARE1 type [11] and 
the other cuts and parameters used are described in the caption of Table 1. The results 
of the various programs agree again within their statistical errors, although the accuracy 
reached in this comparisons is much lower especially at the hard end of the spectrum, due 
to the long time necessary to run OLDBIS to obtain the required accuracy. 

The program BHAGEN-1PH can be used of course also at lower energies, for example 
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in the DA$NE energy range, for which one can simulate the so called large-large angle 
electron-positron double tagging [6]. For the other extremely forward configurations de- 
scribed in [6] the formulae used in BHAGEN-1PH and presented in this paper are not 
adequate. In the Table 2 we show a comparisons of the results of BHAGEN-1PH with the 
results of PHIPHI, which were presented in [6] and [7] . The cuts used are described in the 
caption of Table 2, where the results are seen to agree well, although the large numerical 
errors. Again the accuracy of the comparisons is limited by the accuracy of the PHIPHI 
program results. As the energy is low these comparisons do not test the Z contributions 
to the cross-section and again t-channel contributions are mostly concerned, due to the 
angular cuts used. The tests are completed by the comparison with the BABAMC pro- 
gram [16] for large angles around the Z boson resonance, which allows for the tests of 
s-channel Z boson exchange contributions. As one can see from the Table 3 the agreement 
with BABAMC is good except for the hard part (0.5 Ef, < E 1 < 0.9 Ef,) of the photon 
spectrum at beam energy Eb = 47.5 GeV, where we have found the difference of 0.60(17)% 
. We guess that this difference is due to the approximation used in BABAMC to simulate 
this part of the spectrum, although the accuracy seems to be sufficient for every expected 
experimental precision in this part of the allowed phase space. Finally in Table 4 we 
present some cross-section results obtained with the longitudinally polarized initial elec- 
tron, with beam degree of polarization Pl = ±1, in the experimentally accessible region 
0.1 < E 1 /Eb < 0.5, cuts and parameters are specified in the table caption. 

4. Conclusions. 

We have obtained simple analytic formulae for the cross-section of the e + e~ i— > e + e~7 
process, with Z exchange included and longitudinally polarized initial electron. A fast 
Monte Carlo event generator simulating the process, BHAGEN-1PH, has been imple- 
mented to cover the experimental event selection, although in its present version the 
program is not suitable for configurations where final fermions are extremely forward 
(9± < lmrad). 

Apart from the intrinsic interest for a direct independent test of the electroweak cou- 
plings, the bremsstrahlung process constitutes an important background to other processes. 

Together with the existing soft and collinear corrections the program is used [10,11] to 
calculate radiative corrections to the Bhabha scattering, including Z exchange, for different 
detector arrangements. 
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A 


OLDBIS 


BHAGEN-1PH 


BG/OB 


0.01 - 0.1 


12.235(4) 


12.2374(5) 


1.0002(4) 


0.1 - 0.3 


4.861(4) 


4.8557(3) 


0.9989(9) 


0.3 - 0.5 


1.631(4) 


1.63236(15) 


1.001(3) 


0.5 - 0.7 


0.630(4) 


0.63183(9) 


1.003(7) 


0.7 - 0.9 


0.303(4) 


0.31061(5) 


1.025(14) 



Table 1. Total cross-section (in nanobarns) are compared between OLDBIS [12] 
and BHAGEN-1PH for t-channel contribution to the e + e~ h- > e + e~7 process, for different 
ranges of the photon emitted energy (A = -g 1 ), for the electron and positron angular 
range 3° < 9± < 8° and for beam energy Eb = 47.585 GeV. No cuts are applied on photon 
angles. The statistical variance (one standard deviation) is given in brackets In the third 
column is the ratio of BHAGEN-1PH cross-section to OLDBIS cross-section. 



E™ in (MeV) 


PHIPHI 


BHAGEN-1PH 


BG/PH 


2 


15.9(6) 


16.18(7) 


1.02(4) 


5 


13.0(5) 


13.09(6) 


1.01(5) 


10 


10.7(4) 


10.78(5) 


1.01(4) 


50 


5.66(21) 


5.68(2) 


1.00(4) 


100 


3.6(1) 


3.689(14) 


1.02(3) 



Table 2. Total cross-section (in microbarns) are compared between PHIPHI (values 
in [6] and correction factor in [7]) and BHAGEN-1PH for the e + e~ i— > e + e~7 process, 
for different ranges of the photon emitted energy {E™ in < E 1 < Eb — m 2 e /Eb), for the 
electron and positron angular range 8.5° < 9± < 171.5° and for beam energy Eb = 0.51 
GeV. The minimal allowed angle between photon and final lepton is 0.5°. The statistical 
variance (one standard deviation) is given in brackets. In the third column is the ratio of 
BHAGEN-1PH cross-section to PHIPHI cross-section. 
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E h (GeV) 


A 


BABAMC 


BHAGEN-1PH 


BG/BC 


43.5 


0.1 - 0.5 


74.143(15) 


74.150(5) 


1.0001(3) 


43.5 


0.5 - 0.9 


15.970(13) 


15.9824(14) 


1.0008(9) 


45.5 


0.1 - 0.5 


241.24(8) 


241.29(3) 


1.0002(5) 


45.5 


0.5 - 0.9 


55.34(7) 


55.43(1) 


1.0016(14) 


47.5 


0.1 - 0.5 


63.780(11) 


63.758(5) 


0.99966(25) 


47.5 


0.5 - 0.9 


8.325(13) 


8.3750(8) 


1.0060(17) 



Table 3. Total cross-section (in picobarns) are compared between BABAMC [16] 
and BHAGEN-1PH for the e + e~ ^ e + e~7 process, for different ranges of the photon 
emitted energy (A = -^ L ) and beam energy E b , for the electron and positron angular 
range 40° < 9± < 140°. No cuts are applied on photon angles. The statistical variance 
(one standard deviation) is given in brackets. In the third column is the ratio of BHAGEN- 
1PH cross-section to BABAMC cross-section. The following parameters are used: Z boson 
mass Mz = 91.175 GeV, Z boson width V z = 2.3355 GeV and weak mixing angle sin 6V = 
0.2247. 



E b (GeV) 


A 


Pl = +1 


Pl = 


Pl = -1 


43.5 


0.1 - 0.5 


68.095(3) 


74.150(4) 


80.198(5) 


45.5 


0.1 - 0.5 


200.90(1) 


241.29(2) 


281.69(2) 


47.5 


0.1 - 0.5 


57.469(3) 


63.760(4) 


70.062(5) 



Table 4. Total cross-section (in picobarns) from BHAGEN-1PH for the process 
e + e~ 1— > e + e~7, for longitudinally polarized electron beam with degree of polarization Pl, 
for 0.1 < A < 0.5 (A = ^) and for different beam energies E b , for the electron and 
positron angular range 40° < 9± < 140°. No cuts are applied on photon angles. The 
statistical variance (one standard deviation) is given in brackets. The parameters are as 
in Table 3. 
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